Photon blockade induced by atoms with Rydberg coupling 
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We study the photon blockade of the correlated two-photon scattering in a one-dimensional waveg- 
uide, which couples to two atoms with the Rydberg interaction. We obtain the analytic solution 
to the scattering problem of photonic wave packets on the two-coupled atoms. We examine the 
photon correlation by addressing the two-photon relative wave function and the second-order corre- 
lation function in the single- and two-photon resonance cases. It is found that, in the single-photon 
resonance condition, photon bunching and antibunching exist respectively in the transmitted and 
reflected two photons. In particular, the induced bunching and antibunching become stronger with 
the increase of the inter-atom coupling strength. In addition, we find a phenomenon of bunching- 
antibunching transition caused by the two-photon resonance. 

PACS numbers: 42.50.Ar, 42.50.Pq, 42.79.Gn, 42.65. .k 
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I. INTRODUCTION 

Controllable transport of photons in one dimension 
(ID) waveguides is crucial for realizing all-optical quan- 
tum devices, which are basic elements for implementation 
of photonic quantum information processing. As an ex- 
ample, single-photon transistor [lj, controlling the prop- 
agation of signal photons by another single photon, can 
be used to implement a single-photon quantum switch. 
Up to now, there have been some proposals for coherent 
transport of single- and two-photons in wave guid es with 
linear and nonlinear dispersion relations |2h18|. Most 
of them are based on the physical mechanism of photon 
scattering by tunable targets. 

Two-photon transport in a waveguide is a nontriv- 
ial mission since nonlinear interaction between pho- 
tons fl9l . |20| may be induced by nonlinear scattering tar- 
gets such as a nonlinear cavity [l2j and strongly interact- 
ing atoms 21]. In particular, there exists a kind of non- 
linear interaction could lead to photon blockade [ljl [2(| , 
which can be used to realize a single-photon source. Be- 
cause the photon nonlinearity in a cavity with atoms 
is a second-order effect of the interactions between cav- 
ity field and atoms it is quite natural to ask how 
the atoms could mediate such a nonlinearity for pho- 
ton blockade. Recently, Shi et al. [22j have considered 
two-photon scattering in a waveguide coupled to a cav- 
ity containing a two-level atom. Their results on photon 
blockade can well fit the experiment data [20j |. 

We note that the blockade effect can also occur among 
atoms. In atomic blockade, double excitation of atoms is 
strongly suppressed by some inter-atom coupling [23l - l28j . 
Physically, when the interaction between excited atoms is 
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strong enough, it will be difficult to simultaneously excite 
two atoms due to the coupling induced elevation of total 
energy. Since the simultaneous excitation of two atoms 
needs to absorb two photons, the atomic blockade will 
suppress the two-photon absorption one after another. 
In other words, the atomic blockade can induce the pho- 
ton blockade under certain conditions. With this moti- 
vation, we will investigate in this paper how the atomic 
blockad e l24l . l28j exert an influence on the photon block- 
ade [Hilar 

Deeply understanding the relation between these two 
blockade effects (atomic blockade and photon blockade) 
will provide a straight forward way to probe the nature 
of the inter-atom coupling through photon blockade. It 
is also possible to display the atomic blockade effect due 
to the photon blockade phenomenon. To illustrate the 
physical mechanism behind these scientific and techni- 
cal issues, we study a hybrid system composed by a ID 
linear waveguide, which contains two inter-coupled two- 
level Rydberg atoms. To understand the single photon 
contributions in the two photon process, we first calcu- 
late spatial evolution of the single photon wave packet in 
ID w aveg uide. The similar approach has been used in 
Refs. [lOj, ll2j to study the spatial wave packet scattering 
problem. Then we calculate the evolution of the two- 
photon spatial wave packet in a ID waveguide. Since 
the second-order correlation function can describe 
photon blockade effect, we calculate for both the re- 
flected and transmitted photons respectively. We find the 
second-order correlation function g^ is proportional to 
the joint probability of the photons with some distance. 
We can predict the photon statistical properties by the 
wave function of the two photons. Our calculation shows 
the spatial distribution of the two photons is strongly 
dependent on the inter-atom coupling. The transmit- 
ted photons are strongly bunched except for some spe- 
cial inter-coupling values, which depends on the initial 
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conditions. The reflected photons are strongly repulse 
each other, namely the photon blockade effect, when the 
value of inter-atom coupling is away from the resonance, 
which means the inter-atom coupling equals the sum of 
each photon's center detuning. On the contrary, at this 
resonance, the reflected two photons exhibits bunching 
behavior. We can use this sole bunching behavior for the 
reflected two photons to probe the nature of the inter- 
atom coupling. 

This paper is organized as follows. In Sec. [Ill we intro- 
duce the physical model and its Hamiltonian. In Sees. IIIII 
and IIV1 we consider the single- and two-photon scatter- 
ing with the Laplace transformation, respectively. In 
particular, we describe the two-photon correlation in the 
two-photon relative coordinate space, and calculate the 
second-order correlation function in two-photon reflec- 
tion and transmission. Finally, we draw our conclusion in 
Sec. IVland present the detailed derivation of two-photon 
solution in Appendix. 



II. MODEL SETUP 

We start by considering a ID linear waveguide, which 
contains two two-level atoms coupled via the Rydberg 
interaction (see Fig. [IJ. The model Hamiltonian (with 
h = 1) of the system reads as 



H 



dkLJk{rlr k + llh) + -z-{di + 02) 



-9o 



dk 



1=1,2 



'crt{r k + i k ) + {ri + P k )c 



+£|e)i(e|i (g) |e) 2 <e| 2 . (1) 

The first line of Eq. ([1]) is the free Hamiltonian of the 
fields and atoms. The creation (annihilation) operators 
f]. (r/c) and l\ (l k ) describe, respectively, the right- and 
left-propagating fields in the waveguide, with wave num- 
ber k and frequency Lo k = v p k (hereafter we take light ve- 
locity v p — 1). Pauli's operators a^' y ' z [of = 5 (o^iier^ )] 
are introduced to represent the £th (£ = 1, 2) atom with 
energy level spacing luq. The Hamiltonian in the second- 
line of Eq. ([1]) depicts the atom-field interaction with the 
coupling strength g^. In addition, the last term in Eq. ([1]) 
stands for the Rydberg interaction of strength £ between 
the two atoms [29|, |3(| • 

By introducing the even- and odd-parity modes 

h = ^ k + l k ), c k = -^(fjt - h), (2) 

the Hamiltonian (TTJ can be decomposed into two parts, 
H = + H < e ) with H 

= 



J °° dkw k c\c k and 



dkuj 
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+g / dk[(a+ + a+)b k + b k (a 1 + «r 2 )] 
Jo 

+£|e)i(e|i <8 |e> 2 <e|a, 




FIG. 1. Schematic diagram of the physical setup. Two inter- 
acting atoms are placed in a ID linear waveguide. Photons 
injected from the left-hand side of the waveguide are scattered 
by the atoms. 



where g = \/2go- Obviously, the odd-parity modes de- 
couple with the atoms so that their evolution is free. In 
the following we will mainly deal with the evolution of 
the even-parity modes. 

In a rotating frame with respect to 



H { o e) = 



dkblb k + -£(&i + 6i 



(4) 



the Hamiltonian becomes 

H { j e) = / dkA k b{b k + £|e)i(e|i ® |e) 2 <e| 2 
Jo 
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+g / dk[(a+ + a+% + b^ + a^)], (5) 



where A k ^ uj k — ujq. Based on Hamiltonian ([5]), we can 
work out the photon scattering solution for the system. 
To comprehensively understand the physical picture for 
the photon scattering, we will firstly consider the single- 
photon scattering, and then we will study how does the 
Rydberg interaction affect the two-photon scattering. 



III. SINGLE-PHOTON SCATTERING 

In the following discussion, we employ the time- 
dependent approach to the wave packet scattering [l(| E3 
other than the approach based on stationary state [2j- 
H, S El El to study our photon scattering problem. 
We note that the total excitation number operator N ~ 
d^d + |e)i(e|i + |e) 2 (e| 2 in this system is a conserved 

quantity due to [N, H^] = 0. Consequently, the Hilbert 
space of the system can be divided into the direct-sum 
subspaces with different excitations. For single-photon 
scattering, it is sufficient to consider the scattering prob- 
lems within the single-excitation subspace. In this sub- 
space, an arbitrary state of the system can be expressed 
as 



|*(t))=ai(t)|0)|e) 1 | ff ) a +a a (t)|0)| fl ) 1 |e) 2 
dkp k (t)b{\<b)\g)i\g) 2 , 



(3) 



(0) 



where |0) is the vacuum state and |lfc) represents the 
state with a single photon in the fcth even mode of 
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the waveguide. The variables a\(t), (t), and ft k (t) 
are probability amplitudes. In terms of the Schrodinger 

equation i(d/dt) = flj e) |$(i)), we obtain the fol- 

lowing equation of motion for probability amplitudes 



di(t) 
d 2 (<) 



-ig I dkPk(t), 
lo 



$k(t) = -iA k pk(t) - ig[ai(t) + a 2 (i)] 



(7) 



For single-photon scattering, we assume that the two 
atoms are initially in ground state and a single even- 
mode photon is in a Lorentzian wave packet. The initial 
condition reads 



<*i(0) =0, 
(3 k (0) = Gi 



a 2 (0) = 0, 

„iA fc Z 



A k -S- 



l > 0, 



(8) 



where G\ = Wc/n is the normalization constant, and I is 
the initial distance between the wavefront of the photon 
wave packet and the atoms. In addition, 5 and e are, 
respectively, the center detuning and width of the wave 
packet. The transient solution of these probability ampli- 
tudes might be obtained by the Laplace transformation 
method. For the motivation of scattering problem, we 
focus on the long-time solution, 



Ot\{t 

P k (t 



oo) = 0, a%{t — > 
oo) =i k Pk(0)e~ tAk \ 

where we introduce the phase factor 

A fc 



oo) = 0, 



tk 



A, 



■ vy 
«7 



(9) 



(10) 



with 7 = 2ng 2 being the decay rate of the atom. 

In practice, a single photon should be injected in right- 
or left-propagating modes. Hence, we need to consider 
the single-photon scatting in modes f k and l k . For the 
right-propagating injection, the initial state is 



1-0(0)) 



dk0 k {O)rl\Q) 



dkp k (Q)(b\+cl 



(11) 



V2 Jo 

The wave function in position space reads 

(a#(0)) ^ -iV^G^-^+^ei-x - I), (12) 

where 9{x) is the Heaviside step function. In the long- 
time limit, the single-photon wave function is 



|V>(t->oo)) = / dkl3 k (0)e 
Jo 



(T k 7 



(13) 



where the transmission and reflection amplitudes are ob- 
tained as 

A fc 



T k = 



R k = 



-vy 



(14) 



A fc + 17 ' Afc + ^7 ' 

The output wave function in position space is 

(x\ip(t -> oo)) = ip r (x, I, S) + tpi(x, I, S), (15) 

where 

Vv(z, I, S) = -iV^G^s^e^-^^+^ei-x +*-/), 
ipi{x, 1, 5) = -i^G x R 5 -iee- i{5 - ie){ - x+t - l) B{x + t-l). 

(16) 

We can see from Eq. (I16[) that the waveform (propor- 
tional to \ip r (x, I, S)\ 2 and \ipi(x, 1, 5)\ 2 ) of the transmit- 
ted and reflected photon is the same as that for the in- 
put state (fT2"j) , whereas the amplitudes are normalized by 
the transmission and reflection coefficients (|T,5_i e | 2 and 
\Rs^\ 2 )- 

IV. TWO-PHOTON SCATTERING 

We now turn to study the two-photon scattering. We 
will use the Laplace transformation to obtain the long- 
time state of two photons, and special attention will 
be paid to clarifying the relationship between the atom 
blockade and the photon blockade. 

A. Equations of motion and solution 

In the two-excitation subspace, there are four types 
of basis states: two excitations in atoms (with basis 
1 0) | e) i |e) 2) , one in atoms and the other in fields (with 
bases |lfc)|e)i|<?)2 and |lfc)|g)i|e)2), two in fields (with 
basis \l p , 1 9 ) 1 17)1 1 3)2). Then an arbitrary state in this 
subspace is written as 



| V (t))=A(t)|0)|e) 1 |e) 2 



dfcB fc (t)|lfc)|e)i|5) 2 



dfcC fc (t)|l fc )|0>i|e>2 



dp / dgX> M (i)|lp,lg)|5>ib)a, (17) 
10 Jo 

where A(t), B k (t), C k (t), and D p>q (t) are probability am- 
plitudes. From the Schrodinger equation, we get the fol- 
lowing equations of motion for probability amplitudes, 

/>oo 

A(t) = -i£A(t) - ig / dk[B k {t) + C k (t)}, 
Jo 

poo 

B k (t) = -iA k B k (t) - igA(t) -ig dpD p<k (t), 

Jo 

poo 

C k (t) = -iA k C k (t) - igA(t) -ig d P D p , k (t), 

Jo 

D P , q {t) = -i(A p + A q )D p Jt) - ig[B p (t) + B q (t)} 

-ig[C p (t) + C q (t)}, (18) 
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where we have used the symmetry relation Dk, p (t) = 
D p , k {t). 

We consider the initial state where the two atoms are 
in their ground state and the two photons are in a wave 
packet. The corresponding initial condition reads 



A(0) = 0, 
D M (0) = G 2 



B k (0) = 0, 



a iA p h 



C fc (0) 



0. 



P — Si 



A 



ie A 

e 



q- 02 
>A p l 2 



Si + ie A p — S 2 + ie 



(19) 



where lg is the initial position of the £th {I = 1,2) photon. 
Without loss of generality, we assume 1% > 0, l 2 > 0, and 
l± > l 2 below. The normalization constant reads 



G 2 



4e 2 e -2e(/ 1 -/ 2 ) 

{Si - S 2 f + Ae 2 



-1/2 



(20) 



Under the above initial condition, the long-time so- 
lution for these probability amplitudes might be ob- 
tained with the Laplace transformation, which satisfy 
A(t -> oo) = 0, B k (t -> oo) = 0, C k (t -> oo) = 0, 
and 

D M (t -»■ oo) = (tp*,I>P,«(°) + Jp, g )e- i ( A " +A ")*, (21) 
where i p (t g ) has been defined by Eq. (fl"U|) and 



J P , 9 = 4G 27 5 



i(A p +A,)ii 



(A p + A g - 20 



(A p + i 7 )(A, + 17) (A p + A g - £ + 17) 
1 



(Ap + A, - Si - 5 2 

e -(,i5 2 +e)(h-h) 



2ie) 



-(it 2 +e)(h-l 2 ) 



17 + 02- 

-idi-h) 



it 



- Ly q — 02 



+ ie + 17 



17 - 



+ + 17 



(22) 

The first term in Eq. (|2ip describes two-photon indepen- 
dent scattering process, while the second term represents 
photon correlation induced by scattering. It should be 
pointed out that, when £ = 0, J Pj9 ^ 0. This fact means 
that the photon correlation can be induced even in the 
absence of the Rydberg interaction. 



B. Two-photon wave functions in real space 

We consider a realistic case with the two photons in- 
jected from the left-hand side of the waveguide. Then, 
the initial state of the photons can be written as 



dp / dqD Pi q(0)rlr. 



t*t| 



(23) 
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In terms of the basis wave function 



(xx,x 2 \flfl\Q) =N r . 



,iA p ii iA,i 



e^" X2 x 2 ) (24) 



with Mrr = l/(2-\/27r), the wave function in position 
space of the initial state can be written as 

(xi,x 2 \if)(0)) = -ATT 2 NrrG 2 e [l5l+c){xi+h) e [lS2+e)[x2+h) 
x9(-xi - h)9(-x 2 - l 2 ) + (xi «-» x 2 ). 

(25) 

By introducing the center-of-mass coordinate x c = 
(xi+x 2 )/2, relative coordinate x — xi — x 2 , total momen- 
tum E = Si + S 2 and relative momentum S = (Si — S 2 )/2, 
the wave function ([25]) becomes 



(xi,x 2 \i>(0)) =-4n 2 KrG 2 



e iS 1 (x c +x/2+h) e e(x c +x/2+h) 



y e iS 2 (x c -x/2+l 2 ) e e(x c -x/2+l 2 ) 

X 9(-X c + x/2 - l 2 )6{-x c - x/2 - h) 
+(x o -x)] . (26) 

For the special case of l 2 = h, the wave function reduces 
to 

(xi,x 2 \ip(0)) = -%ir 2 N rr G 2 e^ E+2 ^ +l ^ 

xcos(Sx)9(-x c -h-\x\/2). (27) 

In the long-time limit, the two-photon state can be 
expressed as 

\i){t -> Oo)) = \Ar) + \4>rl) + \i>lr) + W) , (28) 

where 



\iprr) 



oo poo 



Jo 

OO pOO 



Jo 

oo poo 



JO 

OO p oo 



Jo 



dpdqD^e-^+^fplm, 
dpdqD^e-^+^Hlflm, 
dpdqD p \ q e-^ +A ^ilil\<!)), (29) 



ith 



D™„ = 



p,q 

XI 
p,q 



-(TpT q D p , q (0) + J p , q /4), 
^(TpR q D Ptq (0) + Jp, q /4), 
±(R p T q D Ptq (0) + J P JA), 
±(RpR g D p , q (0) + J p>q /4). 



(30) 



Here D^~ q and D p l are, respectively, the two-photon 
transmission and reflection amplitudes, while D p l q and 
Dp t q are the one-photon reflection and one-photon trans- 
mission amplitudes. 



C. Two-photon correlation in position variables 

To characterize the photon correlation, it is conve- 
nient to consider the two-photon transmission or reflec- 
tion cases. For the two-photon transmission, the output 
state of the two right-going photons is 

( Xl ,x 2 \iP rr ) = -An 2 N rr G 2 e l( - E - 2 ^-t) 

Xe K§-iWi+l*) e mi-l*)§ rr ( x } ( 31 ) 

with 

$ rr (x) = T Sl ^ e T S2 - le {e lSx 6(~x c -x/2 + t- h) 
x6(-x c + x/2 + t - l 2 ) + x o -a;] 

Sl - te &2 - lt E-2ie-i + i 1 
Xe i(&-ie+i-y)\ x \ e (_ Xc + t _ - ll).(32) 

When l 2 = h, the output state ([3"Tj) reduces to 

(asi,sa|Vrr> = S7r 2 KrG 2 e^ E+2 ^- t+l ^ 

x9(-x c + t-k-\x\/2)(p rr (x), (33) 

where 



4>rr(x) = T Sl -i e Ts 2 - i<i COs(Sx) 



E-2£-2ie &i{ E_ ie+mw] 



E - £ - 2ze + »7 



(34) 



which satisfies <p rr (—x) — (f> rr (x). In the derivation of 
Eq. (|3"Tj) , we have used the condition 7>e. 

We note that Eq. (f33|) is a product of the center-of- 
mass wave function exp[(iE + 2e)(a; c — t + h)] and the 
relative wave function <p rr (x) in the region defined by 
the step function 9(—x c + t — li — \x\/2). As \4> rr (x)\ 2 
is proportional to the joint probability for two photons 
with a separation x, it could be used to characterize the 
spatial statistics of the two photons. For example, the 
peak and dip feature around |(/> rr (0)| 2 , implies photon 
bunching and antibunching, respectively. In Eq. (|34p. the 
first term of 4> rr (x) describes an independent two-photon 
transmission process, while the second term is a two- 
photon correlation induced by scattering. Physically, the 
present system has two resonant scattering conditions: 
single- and two-photon resonances. When the frequency 
of the a single photon matches the energy separation of 
a single atom, i.e. Si = 5 2 = (or E = and 5 = 0), the 
single photon will resonantly excite the atom. We call 
this as the single-photon resonance condition. On the 
other hand, when the total energy of the two photons 
equals to the energy required to excite the two coupled 
atoms, i.e., 81 + 5 2 = £ (or E = £), the two photons 
can resonantly excite the two coupled atoms, even the 
single-photon process could be off- resonant. This regime 
is called as two-photon resonance. 

In the single-photon resonance case, the independent 
two-photon transmission process will be completely sup- 
pressed, and a pure photon-correlation effect can be seen 
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FIG. 2. (Color online) Plot of |0 rr (a;)| 2 vs the scaled param- 
eters f/7 and yx when (a) E/j = 0, 8 — and (b) E/j = 4, 
5 = 0. (c) and (d) are, respectively, the curves from (a) 
and (b) when the parameter £/7 takes several certain values. 
Here, we consider the near monochromatic limit e/7 = 0.01. 



from \(j) rr (x)\ 2 . In Fig. Efa), we plot \<t> r r{x)\ 2 as a func- 
tion of £/7 and jx. When £ = 0, there is no photon corre- 
lation [dash line in Fig. [2fb)]. This result corresponds to 
the fluorescence-complete-vanishing phenomenon found 
in Ref. [l6j]. For a nonzero £/7, we can see clear ev- 
idence for photon bunching [Fig. [21b)]. In particular, 
with the increasing of £/7, |0, r (O)| 2 increases gradually, 
and saturates when (/7> 1. This means that the pho- 
ton bunching becomes stronger for a larger £/j in the 
single-photon resonance regime. 

In the case with single-photon off resonance (e.g., 
E/j = 4), we plot \4> rr {x)\ 2 vs £/7 and 72; in Figs. [SJc) 
and [D(d) . The curves exhibit photon bunching in most 
region of £/7, but there is also some oscillation pattern 
with respect to 72; due to independent photon transmis- 
sion. However, around the two-photon resonance, i.e., 
£ = E, there is a clear evidence for the photon antibunch- 
ing [Fig. HJd)]. This interesting phenomenon of photon 
statistics is induced by the Rydberg coupling. 

By the similar procedure, using the basis wave function 

(xi,x 2 \llil\0) = M U ( e - iA ^e- iA ^ 2 +xi ^ x 2 )(35) 

with Mu = l/(2V27r), the output state of the two left- 
going photons in the long-time limit t — > 00 is 

{xi,xM = -47r 2 M u G 2 e^-^ h+ ^e^ h -^ 

xe i{E-2ie)(-x c -t)^ n ^ (gg) 
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region of the parameter In addition, it is similar to 
the two-photon transmission in the sense that there also 
exists the photon statistics change induced by the two- 
photon resonance. Around £ = E, we see clear evidence 
for photon bunching [Fig. [3jd)]. We find this bunching 
peak is exactly at £ — E for the reflected photons. This 
change of statistical properties can be used to detect the 
form of the interaction between atoms. 
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FIG. 3. (Color online) Plot of \4>u{x)\ 2 vs £/7 and 7a; when 
(a) E/-y = 0, 8 = and (b) E/-y = 4,8 = 0. (c) and (d) are, 
respectively, the curves from (a) and (b) when the parameter 
f/7 takes several certain values. Here, we take e/7 = 0.01. 



with 

*„(ar) = R Sl -ieR6 2 -ie[e- iSx e(x c + x/2 + t- h) 
x9(x c — x/2 + t - l 2 ) + x <h> —x] 
E-2ie-2£ 

-R5 1 -ieR6 2 -ie E _ 2 . e _^ + . ry 

xe z ^- ie+l ^ xl e(x c + t - \x\/2 - h). (37) 

When li = l 2 , the output state of two reflected photons 
becomes 



x9(x c + t-h-\x\/2)(j)ii(x) 



with 



4>ll(x) = Rs 1 -ieRs 2 -ie [C0s(5x) 

I E - 2ie - 2£ 



2 E - 2ie - £ + ?7 



(38) 



(39) 



In Figs. [3ja) and[3jb), we plot \<j)u(x)\ 2 as a function 
of £/7 and 72; in the single-photon resonance E/j = 
and 5/7 = 0. Similar to the two-photon transmission, 
when £ = 0, there is no photon correlation [dash line in 
Fig.[3]Jb)]. For nonzero £, we can see photon antibunching 
[Fig. Cp)]. With the increasing of £/7, the \(j)u{0)\ 2 de- 
creases gradually, and eventually approaches zero when 

£h » 0. 

On the other hand, |</>n(a;)| 2 in the case with single- 
photon off resonance (-E/7 = 4) is shown in Figs. (3Jc) 
and[3Jd). Though there exists some oscillation with re- 
spect to jx, we can still see photon antibunching in most 



D. Second-order correlation function 

We can also present a quantitative description of the 
statistics of the right- and left-going photons using the 
second-order correlation function g^ 2 \ In particular, we 
will only concern the two-photon reflection and trans- 
mission because the photon statistics in these two cases 
makes sense. From Eq. ([25)). the two-photon states are 
defined by 

\4>sa>) = ^J J dxidx2{x 1 , x 2 \ip ss >)'ipi(x l )ipl, (x 2 )|0) 

(40) 

for s, s' — r,l, where the field operators satisfy the 
bosonic commutation relation [ijj s (xi), ipl( x 2)] = S(x± — 
x 2 ). For state (|40|) . the second-order correlation function 
is 



r,l, (41) 



where 

Gf\xx,T) = (Tpssllplixx^iix-L +T)lJj s (x 1 + t)^ s (xi ) \tp ss ) , 

(42) 



with \%l) aa ) = \ip ss )/ v / (ip ss \ip ss ). Then, by Eqs. g0 
and (|4ip. we get 



(2)/ \ \i>as(xi,Xi + t)| 2 J J dXidx 2 \lp ss (xi,X 2 )\' 2 



2 j dx 2 \ip ss (xi,x 2 )\ 2 J dx 2 \ip ss {x 1 +t,x 2 )| 2 ' 

(43) 



for s = r, I. 

For the two transmitted photons, we introduce a new 
frame of reference as Xi = x[ + 1 and x 2 = x' 2 + t. Cor- 
respondingly, the center-of-mass and relative coordinates 
become x' c = (x' x + x' 2 )/2 and x' = x[ — x' 2 , which have 
the same form as those in the old frame of reference of 
x. In the new frame of reference, the state of the two 
right-going photons is 

(x[,x' 2 \^ rr ) = -87r 2 Af rr G 2 e^ E+2 ^<+^ 

xe{-x' c -h-\x'\/2)<t> rr {x'\ (44) 

which is independent of the time t. In the following, we 
will restrict our calculation in the new frame of reference, 
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and omit the superscript '"" for simplicity. After some 
tedious calculations, we get 



g?\r) = F r (x 1 ,r)\cf >rr (r)\ 2 , 



(45) 



where 



F t (xi,t) = M~H{-Xx - h - t)6(- Xi - h) 

r+oo 

x / e~ 2ex \(l) rr (x)\ 2 dx, (46) 



with 



M r = 4e9(- Xl -h- t)9(-Xi - h)e ie(xi+h+ i ; 



x / dxe-~ ztx \(j) rr {x)\ z / dxe- Mx \4> rr {x)Y 

(47 

Similarly, for the two reflected photons, we also in- 
troduce a new frame of reference as x\ = x' x — t and 
x 2 = x' 2 — t, then the state for the two reflected photons 
becomes 

y.e{x' c -h-\x'\/2)Mx'), (48) 

(2) 

and the correlation function g u (r) reads 



gl 2 \r)=F l (x 1 ,r)\Mr)\' 



where 



Fi{x u t) = M l - 1 6(x 1 -h+ t)6{ Xi - h) 

f + 00 

-2ex\ i /_\|2, 



(49) 



x / e-* x \<t>u(x)\ dx, 



(50) 



with 



Mi = 4e6»(xi -h + t)6(xi - h)e- 4t(xi - ll + ^ 

2cti i _/„,\|2 t / „-2exl \|2, 




y-f -i -io 



tff -10 -10 



-^Irh^J-rV 2 #- 0.5 



x / e^l^CO! / e-" x \M*)\ dx. 

Jli—xi Jli—{xi+T) 

(51) 

In Fig. [H we illustrate the second-order correlation 

(21 (2) 

function g r (r) and g\ (r) as a function of the scaled 
parameters £/7 and 77". Here, the parameters are chosen 
the same as those in Figs. [5] and [3] for comparison. In the 
single-photon resonance, i.e., E/j = and S/j = 0, we 
can see that, when £/ 7 ^ 0, c^ 2) (0) > gf' '(r) [Fig-Ha)] 

(2) (2) 

and <? ; (0) < (? ; (r) [Fig. BJc)], which represent pho- 
ton bunching in transmission and photon antibunching 

(2) 

in reflection, respectively. Interestingly, the g\ '(0) ap- 
proaches zero when £/7 ^> 1> which means that the pho- 
ton blockage phenomenon takes place in the reflection 
[Fig-HIc)]. When £/7 = 0, the two photons are reflected 
completely and independently. At the point £/j = 

#, (2) (t) = 1/2 (for Fock state |2)) and g (2] '(r) does not 
make sense. 



15 



gjt -15 -10 



5/y -15 -10 



FIG. 4. (Color online) The second-order correlation function 
ffr (t) and ffj (t) vs f/7 and 7T when (a,c) E/~f = 0, 5 = 
and (b,d) E/j = 4, 6 = 0. 



On the other hand, Figs.HJb) and|4jd) show that the 
two-photon resonance (£ = E 1 ) will induce the photon 
statistics change between bunching and antibunching. 

(2) f2l 

For the two-photon transmission, we see g r (0) > g r (r) 

(2) (2) 

(bunching) in most region of £/7 and g r (0) < <?r ; (r) 
(antibunching) around the two-photon resonance point 
£, = E [Fig. B|b)]. For comparison, the two-photon re- 

f 2) 

flection result is shown in Fig.[4fd). We can see g\ (0) > 

(2) (2) (2) 

g\ '(t) (bunching) around £ = E and g\ (0) < g\ (r) 

(antibunching) in other region. These results are con- 
sistent with our analysis on photon correlation given in 
Sec Etc). 



V. CONCLUSION 

To conclude, we have studied the transport properties 
of the photonic wave packet, controlled by the Rydberg 
interaction between two atoms in one continuum dimen- 
sion. We find that the quantum statistical properties of 
the scattered photons can be predicted from the relative 
wave function of two photons as well as the second-order 
correlation function. With detailed calculations about 
the relative wave functions, we point out how to control 
the quantum statistical properties of the scattered pho- 
tons in the confined system. To change the quantum sta- 
tistical properties from bunching to antibunching or vice 
versa, for a certain Rydberg interaction strength between 
the two atoms, we adjust the detunings of frequency cen- 
ter of the photonic wave packet from each atom energy 
level spacings, respectively. On the other hand, for a cer- 
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tain detunings of frequency center of the photonic wave 
packet from each atom energy level spacings, we can reg- 
ulate the Rydberg interaction strength by varying the 
distance between the two atoms or select two different en- 
ergy levels of the Rydberg atoms. We can use the change 
of the quantum statistic properties to detect the detail 
fashion of the Rydberg interaction between the Rydberg 
atoms. 
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By eliminating other variables, we obtain the following 
equation for the variable B k (s) as 



[A k -i(s + j)]B k (s) 
1 



2.9 2 
+ig 



dA r 



— OO 

oo 



D P ,k(P) 



A p + A fc 



is 



1 



A p + A k 
-dA v . 



is 



B p (s) 



(A.2) 



Appendix: Solution of Eq. (|18p by the Laplace 
transformation 



In this appendix, we present a detailed derivation of 
Eq. (|2ip by employing the Laplace transformation. Un- 
der the initial condition (fl9|). the equation of motion (fT8|) 
becomes 



The solution of B k (s) can be obtained as 



sA(s) = -itA(s) - iff jf" dk[B k (s) + C k (s)} : B k (s) = Afc ^+ 7 ) [/l + 2e ~ Sh ^ + /a)] 

sB k (s) = -iAMs) - Us) - ,g r d P D pAs) , xkhMhMh - h), (A.3) 

J 

poo 

sC k (s) = -iAfeCfe(s) - igA(s) - ig dpD p , k (s), 



sD M {s) = D p>q (0) - i(A p + A q )D p , q {s) 

-ig[B p (s) + B q (a) + C p (s) + C q (s)]. (A.l) with 
I 

f: 



1 (A fe - 5 2 +ie)(A k +<*i-i(s + e)) (A fe - S x + ie) (A fc + 6 2 - i (s + e)) ' 

e — i(h— (a)A fe g— (ii-i 2 )7 f \ \ 



h 
h 



{A k +ij)(A k -S 2 +ie)(s + e + i(A k + S 1 )) (7 - e - itf 2 ) (s + 7 + e + V-7 + »A fe s + 7 + «f/ 5 

g— (Ji— i2)(*<S 2 +e) ^ 



s + z<5i + ^2 + 2e \ (^7 + 5 2 — ie) (<5 2 - Afe - ie) (s + 7 + i£ 2 + e) (s + ^2 + iAfe + e) 
1 1 



(7 - iS 2 - e) (s + 7 + if) (s + 7 + z<5 2 + e) (s + 7 + if) 



(A.4) 



Then we have The transient solution of D p ^ q (t) can be obtained by the 

inverse Laplace transformation. For studying photon 

r> (r\\ — oi [r ( \ -x- r I w scattering, it is sufficient to get the long-time solution 

D Ptq (s) = PA) r ( TiiX ■ ^ g-en in Eq.©. 

s + t(A p + A q ) 
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